Resonance Scattering in Optical Lattices and Molecules: 
Interband versus Intraband Effects 
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We study the low-energy two-body scattering in optical lattices with higher-band effects included 
in an effective potential, using a renormalization group approach. The approach captures most 
dominating higher band effects as well as all multiple scattering processes in the lowest band. For 
an arbitrary negative free space scattering length(a s ), a resonance of low energy scattering occurs 
as lattice potential depths reaches a critical value v c ; these resonances, with continuously tunable 
positions v c and widths W , can be mainly driven either by intraband or both intra- and interband 
effects depending on the magnitude of a s . We have also studied scattering amplitudes and formation 
of molecules when interband effects are dominating, and discussed an intimate relation between 
molecules for negative a a and repulsively bound states pioneered by Winkler et al.[20]. 
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For a dilute ultracold atomic gas, the two-body s-wave 
scattering length a s is known to be conveniently tun- 
able via magnetic-field- induced Feshbach resonances [1, 
2]. Experimentally in the presence of external trapping 
confinements, however, binary atomic collision proper- 
ties can be dramatically modified as revealed both the- 
oretically and experimentally in three-dimensional (3D) 
harmonic traps [3-5], and in waveguides [6-9]. Remark- 
ably, the waveguide confinement can result in very pe- 
culiar effective potentials as pointed out in a few early 
papers[6, 7]; especially, Olshanii et al. systematically 
studied scattering between atoms in a ID waveguide and 
found that the effective potential for atoms in the lowest 
transverse mode can reach the hardcore limit. Interact- 
ing atoms in optical lattices are another subject that has 
attracted enormous interests for the past few years[10- 
12]. However, till now what happens to binary collisions 
in an optical lattice on the other hand have not been 
thoroughly studied and the subject of molecules of Bloch 
waves is also not well understood. It is becoming essential 
to understand the fundamentals of two-body scattering 
and other few-body physics of Bloch states; such analyses 
should form building blocks for future many-body theo- 
ries and set potential references for quantitative calcula- 
tions of parameters in many-body Hamiltonians. Studies 
of this issue can further cast light on dynamics of collid- 
ing atoms or condensates in optical lattices and coherent 
control of atoms in optical lattices[13]. 

Low energy scattering in an optical lattice was previ- 
ously investigated and resonance scattering was pointed 
out for attractive interactions [14]. Studies there were 
carried out in an approximation where laser potentials 
are approximated as harmonic ones so that the center- 
of-mass motion is decoupled from the relative motion of 
two scattering atoms; effectively the problem was reduced 
to two-body scattering within an individual lattice site, 
which is justifiable for deep lattices. In this Letter to 
reveal how Bloch waves are scattered in optical lattices 



at different depths, we propose an approach to resonance 
scattering without utilizing the approximations of sepa- 
rable potentials in Ref.[14]. Our approach captures most 
dominating higher band effects as well as all intraband 
scattering within the lowest band. It is valid for studies 
of resonances in deep lattices at small a s as well as res- 
onances in shallow ones at large a s . And when applying 
our approach to lower dimensional waveguides, we obtain 
identical results discussed previously [6, 7]. 

Our main proposal is to evaluate U e ft, the effective 
potential for atoms in the lowest band that takes into 
account multiple virtual scattering processes involving 
higher bands, and then apply the same procedure to cal- 
culate the full T-matrix of low-energy scattering. When 
free space scattering lengths \a g \ are comparable to or 
larger than the lattice constant a^, virtual scattering to 
higher bands contributes substantially to scattering in 
the lowest band and resonance scattering is driven by 
both interband and intraband effects. We find that the 
higher-band effects on physical quantities are most pro- 
nounced in shallow lattices near resonances(see Fig. 3, 4). 
When magnitudes of a s are arbitrarily small, resonance 
scattering at the bottom of the lowest band is predom- 
inately driven by intraband virtual scattering and is in- 
duced mainly by the enhanced effective masses of atoms 
in optical lattices. 

To facilitate discussions on low-energy scattering 
atoms, we start with a two-body Hamiltonian 



H = Y J t a \a)(a\+Y J U a p\a)(P\ 1 
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with \ct(/3)) being arbitrary two-body scattering states. 
For scattering in free space with a short range poten- 
tial approximated as U(r) = UoS 3 (r), \a), \[3) = |k, — k) 
and U a p = %■ with f2 the volume; to obtain an effec- 
tive low-energy Hamiltonian, we employ the momentum- 
shell renormalization group(RG) equation approach[15]. 
The key idea here is at an arbitrary cutoff momen- 
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turn A, we can further divide the k-space into two re- 
gions, i.e., a core region defined by |k| < A — 8 A and 
a shell A — SA < |k| < A. Correspondingly, we split 
the Hamiltonian at a given cutoff A into three pieces 
-ff(A) = H < + H > + H >< which respectively describe 
interacting atoms within the core, within the shell and 
the scattering in between. For atoms with |k| <C A, the 
second-order virtual scattering into high energy states 
within the shell caused by H >< modifies the low-energy 
scattering amplitudes and results in a correction (5U) in 
H < . One can then obtain a differential RG equation for 
effective potential U(A) in if (A) in terms of the cutoff A, 
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So the effective potential U for scattering atoms at mo- 
menta smaller than A is renormalized due to the coupling 
to virtual states at larger momenta and is given as 
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U(A) U(A*) 



Q ^ 2e k 

A<|k|<A* 



(3) 



Boundary conditions U(A*) = Uq and U(0) = To relate 
Uq to the low-energy scattering length a s (= ^r 2 -) via 

i^k = m + n£|k|<A* ab- (A * is an ultraviolet m °- 

mentum cut-off that is set by the range of interactions.) 
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FIG. 1: (color online). Virtual scattering processes that con- 
tribute to scattering potential To in optical lattices. Differ- 
ent two-particle vertex functions are indicated below the di- 
agrams. Solid (dashed) lines between two vertices represent 
states in the lowest n — band (higher n > bands). 

One can apply a similar idea to optical lattices. 
For a cubic optical lattice with potential V(r) = 
V J2i= x y z sin 2 (^p) and spacing a L — -|, Bloch states 
</>nk(r) = -j= J2g w "( k + G)e j(k+G) ' r and energies e nk 
can be numerically obtained; here n are band indices; 
k G £>Z(Brillouin zone) is a quasi- momentum and G 
(and Q below) are the reciprocal lattice vectors; eigen- 
vector w„(q) is the Wannier wavefunction in momentum 
space. We specify potential depth Vb in units of the recoil 



energy E R 



via a dimensionless quantity v 



First we calculate the interaction matrix element be- 
tween Bloch states \a) = |{m, k}, {n, — k}) and |/3) = 



|{m',k'},{n',-k'}>, 

U af} ee (a\U\(3) =^J2 M « *K> ( 4 ) 
Q 

here = £ G w m(k + G)w n (Q - k G). Relevant 
matrix elements of U a p can be classified into three cat- 
egories: A) |{0, k}; {0, -k}) O |{0, k'}; {0, -k'}), i.e. 
scattering within the lowest band, U a p are given as 
M[/o[16]; these represent the most dominating processes; 
B) \{n, k}; {n, -k}) <-> \{n', k'}; {n', -k'}) with n ^ or 
n' 7^ which constitute the most important scattering 
processes involving higher bands, give the next dominat- 
ing contributions that are approximately equal to Uq (de- 

2 

viations are typically of order of ^ or l ess m shallow lat- 
tices); C) |{n,k};{n,-k}) o |{m', k'}; {n 1 , -k'}) with 
m! 7^ n', i.e. scattering involving two atoms in different 
bands; they contribute the least in shallow lattices (of or- 
der of | or less) because of the approximate translational 
symmetry. 

In shallow lattices at an arbitrary a s , we can always 
neglect matrix elements in C-class and only keep those in 
A- and B-class. In deep lattices near resonances where a s 
are small, we keep matrix elements in B-class to remove 
the ultraviolet divergence when summing up the virtual 
scattering to high energy states; the residue higher band 
effects after regularization turn out to be negligible and 
we again neglect C-class scattering processes; and our 
treatments of scattering processes within the lowest band 
become exact in this limit. However, for large a s and 
deep lattices that are away from the resonances, con- 
tributions from C-class scattering can be comparable to 
other classes; and by neglecting C-class contributions, we 
obtain in this limit estimates only good for qualitative 
understanding. To study resonances, below we adopt a 
simplest two-coupling-constant modcl(See Fig.l) which 
yields reasonable estimates of higher band effects. 

Using the general features of U a p discussed above and 
following the idea outlined before Eq.(3), we obtain the 
effective potential U e ff for the lowest band and further 
calculate the scattering potential To for states near e„ k = 
0, as diagrammatically shown in Fig.l, to be 



1 



mrj 



Ueff 
1 



mrj 



47tolM a 



Co 



(5) 



with C\ 2 defined as 



Ci = 
C 2 = 
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Here r, = (1 + (1 - £„>o,k ^J" 1 is close to 

unity in the regions that interest us[17]; evidently C\ and 



3 



Ci are respectively ascribed to interband and intraband 
scattering effects. Note that when a s is much bigger than 
a L, U e ff saturates at a value of — 47ra L M/mC 1 . 
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FIG. 2: (a)Interband(Ci) and intraband(C2) effects vs lattice 
potentials v. Inset shows the ratio. (b)Resonance position 
v c vs ai/a 3 ; dashed line is v c calculated without interband 
contributions (i.e. G\ = 0.) (c)Width W(v c ). Solid and 
dashed lines are fit to ~ and -^-\pui respectively. 



Our results of C\ t % are shown in Fig. 2. At v — and 
M = 1, we reproduce the free space result Tq — 4ira s /m. 
With increasing v, the intraband scattering gradually 
takes a dominating role over other ones, reflected by a 
much more rapid increase of C2 than C±. For instance 
at v = 5, C1/C2 = 0.21. In the large-w limit, with the 
lowest band spectrum = £$^(1 — cosfcjOj,), t being 
the hopping amplitude, we find that 
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To obtain Bloch wave scattering length auoch we first 
introduce an effective (band) mass m e ff = 1/ g^ fc |o and 
relate it to the scattering potential Tq = ATtauoch/meff- 
For a negative a s , a resonance {auoch — > 00) occurs at 
lattice potential v c when — = — (C2 — Ci). Across the 
resonance, auoch obeys an asymptotic equation 

atioch _ W(v c ) 



(8) 



In the limit of \a s \ -C az,, v c and VF can be estimated 
using Eq.(7); in the opposite limit, they can be obtained 
using the perturbation theory with respect to v, 




\a s \ < a L 
\a„\ > a L 



Both v c and W are continuously tunable by varying a s . 
For ultracold isotopes with negative zero-field scattering 
lengths such as 85 i?6|2,2) (-390a ), 39 ^|1,1) (-45a ) 



and 7 Li\2,2) (— 27ao), using parameters in [18] we find 
resonances at v c — 5.5, 10.0, 11.7 respectively; for 87 Rb 
and 40 K atoms with interspecies scattering length abf — 
— 177ao, resonance scattering occurs at v c = 7.1. 

For very small \a s \, \a s \ <C the effective potential 
U e ff can be simply related to the on-site interaction Uh 
in the Hubbard model, = ^ (O = iV L a|). Follow- 
ing Eq.(5-8), we express auoch as 
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(9) 



which predicts a resonance at jj— s» —0.25. 
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FIG. 3: (color online), (a) Real(/ij) and imaginary(//) com- 
ponents of f{E) = ^^(T- 1 (E) - Tq 1 ) in an optical lattice 
with depth v = 2.5 (the band width 6t = 1.55-E_r). fx is re- 
lated to the normalized density of state p(e), // = 4p(j^-). 
/(_E) calculated without interband effects (Ci = 0) or us- 
ing a Hubbard (Hub) model are also shown. (b,c) Scattering 
cross sections(47ra|,|x| 2 , near the bottom of the band) and 
phase shifts(S) at a 3 = — 0.5dL (resonance position v c = 2.37) 

and v = 2.5, calculated using T-matrix, |xl e * e — "~^£r~~" • 
Dashed lines are the data without interband effects. 

We now turn to scattering matrix T{E) for two atoms 
with total energy £'[19]. Using identical diagrams as 
shown in Fig.l, one can introduce a Lippmann-Schwinger 
equation for two scattering atoms in optical lattices; the 
solutions for T(£) can be obtained as 



£(£) 



mrj(E) a L 
AitcilM a s 



d(E) + C 2 (E)), (10) 



where Ci^(E) can be obtained by substituting e n k in 
of Eq.(5),(6) with e„ k - £/2 - i0+. At small £, the first 
two terms in the bracket in Eq.(10) can be approximated 
to be U e ff that dictates the low energy scattering. Fig. 3 
shows cross sections and phase shifts in shallow lattices 
where higher band effects are dominating (See (b),(c)). 
When £ is approaching zero, asymptotically we have 
T -i {E) = !2gfi(o^+/3fe|o£ + »A: B ),_fc B = ^~E 
and /3 approaches 0.037T when Ci is negligible. T-matrix 
and scattering phase shifts in optical lattices exhibit 
much richer £-dependence than in free space; this is 
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mainly due to a relatively large range of effective inter- 
actions (of order of oll) in the lowest band compared to 
that of free space resonances, or a small resonance energy 
width (of order of t as suggested in Eq.(9)). 




FIG. 4: (color online). Real momentum distribution n^{q x = 
q y = 0, normalized) of bound states at different potential 
depths v(> v c )~, a, s — — 0.5a,L and the resonance occurs at 
v c = 2.37. Right inset shows results without higher band 
effects(i.e. Ci = 0) and the left one shows the binding energy 
(dashed line is the estimate without Ci); note that neglecting 
higher band effects severely overestimates \Eb\ leading to a 
much less singular momentum distribution function n q . 



bound states near the bottom of lowest band for a neg- 
ative a s therefore imply resonance scattering and bound 
states near the the top of the band for a positive scat- 
tering length — a s . Note that our equation for the repul- 
sively bound states in this particular limit is identical to 
the one in Rcf . [20] . 

In conclusion, we have developed an approach to 
low-energy resonance scattering in optical lattices tak- 
ing into account not only the intraband physics but 
more importantly higher band effects. The resonance 
scattering in optical lattices offers an alternative path 
to unitary cold Bose gases so far mainly studied via 
Feshbach resonances [22]. Resonances can also be uti- 
lized to study exciting few-body physics of heteronuclear 
molecules[23] and Efimov states. We thank Immanuel 
Bloch, Hanspeter Biichler, Gora Shlyapnikov, Victor Gu- 
rarie, Maxim Olshanii, Dmitry Petrov, Leo Radzihovsky 
and Ruquan Wang for stimulating discussions and the 
KITPc 2009 cold atom workshop in Beijing for its hos- 
pitality. This work is in part supported by NSFC, 973- 
Project (China), and by NSERC (Canada), Canadian In- 
stitute for Advanced Research. 



Beyond v c , a stable molecule can be formed with a 
binding energy Eb{< 0). Eb can be obtained by solving 
the following two-body equation 



= 



mrj 



(— -C 1 {E B ) + C 2 {E B )) 1 (11) 



Near resonances, \E B | is proportional to a b f och (see Eq.(8) 
for auoch) with the same scaling dimension as in free 
space. For a bound state \^) = Y^ n k c nk.ipnk.^ n -k > n i = 
SnkG <5q,k+G \w n (k + G)| 2 |c„k| 2 and c„k is proportional 
to I/(2e„k — Eb)- In Fig. 4, we plot n q for a s — —O.baL 
where higher band effects are dominating. 

In the limit of deep lattices[20, 21], one can neglect 
higher band effects by setting C X {E) = in Eq.(10),(ff) 
and rj(E) = 1. Both the T- matrix and the binding en- 
ergy in this limit exhibit a generalized particle-hole sym- 
metry due to a property of the single particle density 
of states, p(e) = p(6t — e). So for a given scattering 
length a s , one finds that ReT- % (E) + ReT~ 1 ( 12t - E) = 
i? e r- 1 (I2t) + i?eT- 1 (0) and ImT- l {E) = ImT- 1 (l2t- 
E)(see Fig. 3a). More important, the stable molecules 
below the lowest band for negative scattering lengths 
a s (< 0) have close connections to mid-gap repulsively 
bound states for positive scattering lengths that were first 
thoughtfully pointed out by Winkler et al.[20}. In addi- 
tion, the T-matrix for negative a s can also be related to 
that for positive a s via a simple reflection symmetry. In- 
deed, by examining Eq.(f 0), (11) in the limit of deep lat- 
tices we verify the following exact relations between a s (< 
0) and -a s (> 0) cases, T(E,a s ) = -T*(12t - E, -a s ), 
—Eb{ci s ) = Eb(—ci s ) — 12t; resonance scattering and 
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